arXiv:1503.06214v2 [astro-ph.CO] 18 May 2015 


Emulating the CFHTLenS Weak Lensing data: Cosmological Constraints from 

moments and Minkowski functionals 

Andrea Petri, ^’^’0 Jia Liu,^ Zoltan Haiman,^ Morgan May,^ Lam Hui,^ and Jan M. Kratochvid 

^Department of Physics, Columbia University, New York, NY 10027, USA 
^Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 
^Department of Astronomy, Columbia University, New York, NY 10027, USA 
‘^Astrophysics and Cosmology Research Unit, University of KwaZulu-Natal, Westville, Durban 4000, South Africa 

(Dated: May 19, 2015) 

Weak gravitational lensing is a powerful cosmological probe, with non-Gaussian features poten¬ 
tially containing the majority of the information. We examine constraints on the parameter triplet 
{Dm, w, as) from non-Gaussian features of the weak lensing convergence field, including a set of mo¬ 
ments (up to P** order) and Minkowski functionals, using publicly available data from the 154 deg^ 
CFHTLenS survey. We utilize a suite of ray-tracing N-body simulations spanning 91 points in 
(Dm, w, as) parameter space, replicating the galaxy sky positions, redshifts and shape noise in the 
CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a 
function of (Dm, w, ag), and use it to compute the likelihood function and parameter constraints. 

We employ a principal component analysis to reduce dimensionality and to help stabilize the con¬ 
straints with respect to the number of bins used to construct each statistic. Using the full set of 
statistics, we find Es = as(Dm/0.27)°'^^ = 0.75 ± 0.04 (68% C.L.), in agreement with previous 
values. We find that constraints on the (Dm,as) doublet from the Minkowski functionals suffer 
a strong bias. However, high-order moments break the (Dm,as) degeneracy and provide a tight 
constraint on these parameters with no apparent bias. The main contribution comes from quartic 
moments of derivatives. 

PACS numbers: 98.80.-k, 95.36.-l-x, 95.30.Sf, 98.62.Sb 
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I. INTRODUCTION 

Weak gravitational lensing (hereafter WL) is emerging 
as a promising technique to constrain cosmology. Tech¬ 
niques have been developed to construct cosmic shear 
fields with shape measurements in large galaxy cata¬ 
logues. Although the shear two-point function (2PCF) 
is the most widely studied cosmological probe (see, e.g. 
[H), alternative statistics have been shown to increase 
the amount of cosmological information one can extract 
from weak lensing fields. Among these, high-order mo¬ 
ments (@-0), three-point functions (Q-Q), bispectra 
([iMl), peak counts ((Hill) and Minkowski Func¬ 
tionals (1^,1^) have been shown to improve cosmolog¬ 
ical constraints in weak lensing analyses. 

In this work, we use the publicly available CFHTLenS 
data, consisting of a catalog of ^4.2 million galax¬ 
ies, combined with a suite of ray-tracing simulations in 
91 different cosmological models to derive constraints 
on the cosmological parameters Dm, erg and the dark 
energy (DE) equation of state w. The statistics we 
consider in this work are the Minkowski functionals 
(MFs) and the low-order moments (LM) of the con¬ 
vergence held. Cosmological parameter inferences from 
CFHTLenS have been obtained using the 2PCF [I|, and 
a number of authors have investigated the constraining 
power of CFHTLenS using statistics that go beyond the 
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usual quadratic ones. Fu et al. used three-point cor¬ 
relations as an additional probe for cosmology, and found 
modest (10 — 20%) improvements over the 2PCF. These 
results rely on the third order statistics systematic tests 
performed by Simon et al. [2^. 

Liu et al. [2^ have found a more signiheant (50 — 60%) 
tightening of the D and erg constraints, utilizing the abun¬ 
dance of WL peaks. Cosmological constraints using WL 
peaks in CFHTLenS Stripe82 data have also been inves¬ 
tigated by Liu et. al. (2^. Finally, closest to the present 
paper, Shirasaki & Yoshida [2^ investigated constraints 
from Minkowski Functionals, including systematic errors. 
Our study represents two major improvements over pre¬ 
vious work. First, constraints from the MFs in ref. [2^ 
were obtained through the Fisher matrix formalism, as¬ 
suming linear dependence on cosmological parameters. 
Our study utilizes a suite of simulations sampling the cos¬ 
mological parameter space, mapping out the non-linear 
parameter-dependence of each descriptor. Second, we in¬ 
clude the LMs as a set of new descriptors; these yield the 
tightest and least biased constraints. 

This paper is organized as follows: we first give an 
overview of the CFHTLenS catalogs, and summarize the 
adopted data reduction techniques. Next, we give a de¬ 
scription of our simulation pipeline, including the ray¬ 
tracing algorithm, and the procedure used to sample the 
parameter space. We call the statistical weak lensing 
observables - the power spectrum, Minkowski function¬ 
als, and moments - ’’descriptors” throughout the paper. 
We discuss the calculation of the descriptors, including 
dimensional reduction using a principal component anal- 
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ysis, and the statistical inference framework we used. We 
then describe our main results, i.e. the cosmological pa¬ 
rameter constraints. To conclude, we then discuss our 
findings and comment on possible future extensions of 
this analysis. 


II. DATA AND SIMULATIONS 
A. CFHTLenS data reduction 

In this section, we briefly summarize our treatment of 
the public CFHTLenS data. For a more in-depth de¬ 
scription of our data reduction procedure, we refer the 
reader to (^ . 

The CFHTLenS survey covers four sky patches of 64, 
23, 44 and 23 deg^ area, for a total of 154 deg^. The 
publicly released data consist of a galaxy catalog created 
using SExtractor [2^ , and includes photometric redshifts 
estimated with a Bayesian photometric redshift code [1^ 
and galaxy shape measurements using lensfit [3l|, . 

We apply the following cuts to the galaxy catalog: 
mask < 1 (see Table B2 in [1^), redshift 0.2 < z < 1.3 
(see dH), fitclass = 0 (which requires the object to 
be a galaxy) and weight w > 0 (with larger w indicat¬ 
ing smaller shear measurement uncertainty). Applying 
these cuts leaves us 4.2x10® galaxies, 124.7 deg^ sky cov¬ 
erage, and average galaxy density Ugai ~ 9.3arcmin“^. 
The catalog is further reduced by ~ 25% when one re¬ 
jects fields with non-negligible star-galaxy correlations. 
These spurious correlations are likely due to imperfect 
PSF removal, and do not contain cosmological signal. 
These cuts are consistent with the ones adopted by the 
CFHTLenS collaboration (see [^1. 

The CFHTLenS galaxy catalog provides us with the 
sky position 0, redshift z{6) and ellipticity e{0) of each 
galaxy, as well as the individual weight factors w{0) 
and additive and multiplicative ellipticity corrections 
c{6),m{0). Because the CFHTLenS fields are irregularly 
shaped, we first divide them into 13 squares (subfields) 
to match the shape and Ri 12 deg^ size of our simulated 
maps (see below). These square-shaped subfield maps 
are pixelized according to a Gaussian gridding procedure 

Ef=lWi\0-0MO^)[^ + mi0)] ’ 

where the smoothing scale 6g has been fixed at 1.0 arcmin 
(but varied occasionally to 1.8 and 3.5 arcmin for specific 
tests described below) and m, c refer to the multiplicative 
and additive corrections of the galaxies in the catalog. 

Using the ellipticity grid e{0) as an estimator for the 
cosmic shear we perform a non-local Kaiser- 

Squires inversion [33| to recover the convergence k{0) 


from the E-mode of the shear field. 
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The simulated k maps we create below are 12 deg^ in 
size and have a resolution of 512 x 512 pixels. The 
CFHTLenS catalogs contain masked regions (which in¬ 
clude the rejected fields and the regions around bright 
stars). We first create gridded versions of the observed 
K maps matching the size and pixel resolution of our 
simulated maps, with each pixel containing the number 
of galaxies (jigai) falling within its window. We then 
smooth this galaxy surface density map with the same 
Gaussian window function as equation © and remove 
regions where Ugai < 5 arcmin^ (see [1^). Regions with 
low galaxy number density can induce large errors in the 
cosmological parameter inferences. 


B. Simulation design 

We next give a description of our method to sam¬ 
ple the parameter space with a suite of N-body simula¬ 
tions. We wish to investigate the non-linear dependence 
of the descriptors (in this work, Minkowski Functionals 
and moments of the k field) on the parameter triplet 
p = (flm,re,erg), while keeping the other relevant pa¬ 
rameters (/i, Hb,ns) fixed to the values (0.7, 0.046, 0.96) 
(see M)- We sampled the E-dimensional (E = 3 in this 
case) parameter space using an irregularly spaced grid. 
The grid was designed with a method similar to that used 
to construct an emulator for the matter power spectrum 
in the Coyote simulation suite [s^. Given fixed available 
computing resources, the irregular grid design is more ef¬ 
ficient than a parameter grid with regular spacings: to 
achieve the same average spacing between models in the 
latter approach would require a prohibitively large num¬ 
ber of simulations. 

We limit the parameter sampling to a box whose sides 
range over flm G [0.07,1], w G [—3.0,0], as G [0.1,1.5]. 
These are large ranges, with most of the corresponding 
3D parameter volume ruled out by other cosmological ex¬ 
periments. However, our focus in this work is to quantify 
the constraints from CFHTLenS alone, which, by itself 
has strong parameter degeneracies. We next map this 
sampling box H into a hypercube of unit side. We want 
to construct an irregularly spaced grid consisting of N 
points X,; G [0,1]^. Let a design V be the set of this 
irregularly spaced N points. Our goal is to find an opti¬ 
mal design, in which the points are spread as uniformly 
as possible inside the box. Following ref. [s^, we choose 
our optimal design as the minimum of the cost function 
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This problem is mathematically equivalent to the min¬ 
imization of the Coulomb potential energy of N unit 
charges in a unit box, which corresponds to spreading 
the charges as evenly as possible. Finding the optimal 
design 2 ?^ that minimizes o can be computationally 
very demanding, and hence we decided to use a simplified 
approach. Although approximate, the following iterative 
procedure gives satisfactory accuracy for our purposes: 

1. We start from the diagonal design Vq: xf = i/{N— 
1 ) for d = 

2. We shuffle the coordinates of the par¬ 
ticles in each dimension independently 

xf = T’d ..., l), where 'Pi,...,'Pd are 

random independent permutations of (1, 2,..., N). 

3. We pick a random particle pair {i,j) and a random 

coordinate d £ D} and swap xf Xj. 

4. We compute the new cost function. If the value 
is less than in the previous step, we keep the ex¬ 
change, otherwise we revert the coordinate swap. 

5. We repeat steps 3 and 4 until the relative cost func¬ 
tion change is less than a chosen accuracy parame¬ 
ter e. 

We have found that for = 91 grid points 0(10®) 
iterations are sufficient to reach an accuracy of e ~ 10 “^. 
Once the optimal design 7)^ bas been determined, we 
invert the mapping 11 —>■ [ 0 , 1 ]® to arrive at our simulation 
parameter sampling pg. We show the final list of grid 
points in Table Hand Figured] 

For each parameter point on the grid pg we then run 
an A^-body simulation and perform ray tracing, as de¬ 
scribed in § III C[ to simulate CFHTLenS shear catalogs. 
Throughout the rest of this paper, we refer to this set 
of simulations as CFHTemul. Additionally, we have run 
50 independent A^-body simulations with a fiducial pa¬ 
rameter choice po = (0.26,—1.0, 0 . 8 ), for the purpose of 
accurately measuring the covariance matrices, needed for 
the parameter inferences in 9III Bl This additional suite 
of 50 simulations will be referred to as CFHTcov. 


C. Ray—Tracing Simulations 

The goal of this section is to outline our simulation 
pipeline. The fluctuations in the matter density field 
between a source at redshift z and an observer located 
on Earth will cause small deflections in the trajectories 
of light rays traveling from the source to the observer. 
We estimate the dark matter gravitational potential run¬ 
ning A^-body simulations with N = 512® particles, using 
the public code Gadget2 [s^. We adopted a comoving 
box size of 240h“®Mpc, corresponding to a mass resolu¬ 
tion of 7.4 X lO®h“®M 0 . The simulations include dark 
matter only, and the initial conditions were generated 


with N-GenIC ai z = 100, based on the linear mat¬ 
ter power spectrum created with the Einstein-Boltzmann 
code GAME [s^ . Data cubes were output at redshift in¬ 
tervals corresponding to 80/i“® (comoving) Mpc. 

Using a procedure similar to refs. the equation 

that governs the light ray deflections can be written in 
the form 

^^ = -_V.,4>(x^(x),x), (5) 

where x is the radial comoving distance, xj_ = xP refers 
to two transverse coordinates (with the angular sky co¬ 
ordinates, using the flat sky approximation)a, x(x) is the 
trajectory of a single light ray, and $ is the gravitational 
potential. 

Suppose that a light ray reaches the observer at an 
angular position 6 on the sky: we want to know where 
this light ray originated, knowing it comes from a red¬ 
shift Zs- To answer this question we need to integrate 
equation ([5|) with the initial condition = 0 up to 

a distance Xs = xi^s) to obtain the source angular po¬ 
sition fi{xs\0). Since light rays travel undeflected from 
the observer to the first lens plane, the derivative initial 
condition in the Cauchy problem (O reads ${0\9) = 0. 
We indicate the derivative of with respect to % as fi. 
We use our proprietary implementation Inspector Gad¬ 
get (see, e.g., [^) to solve for the light ray trajecto¬ 
ries based on a discretized version of equation (151) t hat 
is based on the multi-lens-plane algorithm (see [39| for 
example). Applying random periodical shifts and rota¬ 
tions to the A^-body simulation data cubes, we generate 
R = 1000 pseudo-independent realizations of the lens 
plane system used to solve ®. Once we obtain the light 
ray trajectories, we infer the relevant weak leasing quan¬ 
tities by taking angular derivatives of the ray deflections 
A{xs',9) = dfi{xs:9)/d6 and performing the usual spin 
decomposition to infer the convergence k and the shear 
components ( 7 ^, 7 ^), 

A(Xs;9) = (l-K(xs;9))I-x\xs;9)cr^-7^(xs;9)a^ ( 6 ) 

where / is the 2 x 2 identity and ct®’® are the first and third 
Pauli matrices. We perform this procedure for each of the 
R realizations of the lens planes and we obtain R pseudo¬ 
independent realizations of the 7 weak leasing field. We 
use these different random realizations to estimate the 
means (from the CFHTemul simulations) and covariance 
matrices (from the CFHTcov simulations) of our descrip¬ 
tors. Since the random box rotations and translations 
that make up the CFHTcov simulations are based on 50 
independent A^-body runs, we believe the covariance ma¬ 
trices measured from this set to be more accurate than 
the ones measured from the CFHTemul set. 

The convergence k is related to the magnification, 
while the two components of the complex shear 7 = 
7 ® + 17 ^ are related to the apparent ellipticity of the 
source. Given a source with intrinsic complex ellipticity 
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w 

0-8 

1 

0.136 

-2.484 

1.034 

26 

0.380 

-2.424 

0.199 

51 

0.615 

-1.668 

0.185 

76 

0.849 

-0.183 

0.821 

2 

0.145 

-2.211 

1.303 

27 

0.389 

-0.939 

0.454 

52 

0.624 

-2.757 

0.327 

77 

0.859 

-1.182 

1.415 

3 

0.155 

-0.393 

0.652 

28 

0.399 

-1.938 

1.500 

53 

0.634 

-1.575 

0.976 

78 

0.869 

-2.031 

0.227 

4 

0.164 

-2.181 

0.313 

29 

0.409 

-2.940 

0.737 

54 

0.643 

-2.454 

1.444 

79 

0.878 

-2.697 

0.524 

5 

0.173 

-0.423 

1.231 

30 

0.418 

-1.758 

0.383 

55 

0.652 

-1.029 

1.458 

80 

0.887 

-0.363 

0.439 

6 

0.183 

-0.909 

0.269 

31 

0.427 

-2.910 

0.411 

56 

0.661 

-0.486 

0.892 

81 

0.897 

-0.999 

0.468 

7 

0.192 

-1.605 

1.401 

32 

0.436 

-0.060 

0.878 

57 

0.671 

-2.364 

0.793 

82 

0.906 

-1.698 

1.273 

8 

0.201 

-2.787 

0.807 

33 

0.446 

-1.212 

1.486 

58 

0.681 

-2.970 

0.610 

83 

0.915 

-2.544 

1.175 

9 

0.211 

-0.333 

0.341 

34 

0.455 

-2.637 

1.373 

59 

0.690 

-1.332 

0.482 

84 

0.925 

-0.636 

1.259 

10 

0.221 

-1.485 

0.666 

35 

0.464 

-2.121 

0.906 

60 

0.700 

-0.273 

0.283 

85 

0.943 

-2.394 

0.835 

11 

0.239 

-1.848 

0.962 

36 

0.474 

-1.302 

0.114 

61 

0.709 

-2.061 

0.425 

86 

0.953 

-1.545 

0.355 

12 

0.249 

-2.727 

0.369 

37 

0.483 

-1.515 

0.680 

62 

0.718 

-1.728 

1.472 

87 

0.963 

-2.151 

0.510 

13 

0.258 

-1.395 

0.241 

38 

0.493 

-0.243 

0.297 

63 

0.728 

-0.120 

0.596 

88 

0.972 

-0.666 

0.694 

14 

0.267 

-2.667 

1.317 

39 

0.502 

-1.152 

1.189 

64 

0.737 

-2.847 

1.203 

89 

0.981 

-1.242 

1.048 

15 

0.276 

-0.849 

1.429 

40 

0.512 

-0.819 

0.849 

65 

0.746 

-0.090 

1.118 

90 

0.991 

-1.908 

1.020 

16 

0.286 

-1.272 

1.104 

41 

0.521 

-2.334 

0.538 

66 

0.755 

-0.456 

1.359 

91 

1.000 

-1.425 

0.708 

17 

0.295 

-1.878 

0.100 

42 

0.530 

0.000 

0.624 

67 

0.765 

-2.091 

1.076 

- 

- 

- 

- 

18 

0.305 

-0.879 

0.765 

43 

0.540 

-0.030 

1.161 

68 

0.775 

-1.122 

1.132 

- 

- 

- 

- 

19 

0.315 

-2.241 

0.638 

44 

0.549 

-1.818 

1.287 

69 

0.784 

-1.062 

0.779 

- 

- 

- 

- 

20 

0.324 

-2.001 

1.217 

45 

0.558 

-2.577 

1.146 

70 

0.794 

-1.365 

0.156 

- 

- 

- 

- 

21 

0.333 

-0.213 

0.552 

46 

0.568 

-0.516 

1.331 

71 

0.803 

-2.607 

0.255 

- 

- 

- 

- 

22 

0.342 

-2.817 

1.062 

47 

0.577 

-3.000 

0.948 

72 

0.812 

-1.788 

0.722 

- 

- 


- 

23 

0.352 

-0.576 

1.090 

48 

0.587 

-2.304 

0.128 

73 

0.821 

-2.880 

0.863 

- 

- 

- 

- 

24 

0.361 

-0.606 

0.171 

49 

0.596 

-0.696 

0.496 

74 

0.831 

-0.759 

0.213 

- 

- 

- 

- 

25 

0.370 

-0.303 

1.345 

50 

0.606 

-0.789 

0.142 

75 

0.840 

-2.274 

1.387 

- 

- 

- 

- 


TABLE I. List of the CFHTemul grid points in the 3D cosmological parameter space. 
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FIG. 1. and {flm, os) projections of the final simulation design. The blue points correspond to the CFHTemul simulation 

set, which consists of one Ai-body simulation per point, while the red point corresponds to the CFHTcov simulation set, which 
is based on 50 independent A^-body simulations. 
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Gg = Cg + iel, its observed ellipticity will be modified to 


ea+K 

l+g*es 


1+ge* 

ej+g. 


|g| < 1 
|g| > 1 


(7) 


where g = 7/(1 — k) is the reduced shear. 

For each simulated galaxy, we assign an intrinsic ellip¬ 
ticity by rotating the observed ellipticity for that galaxy 
by a random angle on the sky, while conserving its magni¬ 
tude |e|. To be consistent with the CFHTLenS analysis, 
we adopt the weak leasing limit (I 7 I <C 1 ,k <C 1 ), i.e. 
g ft! 7 and e ft -I- 7 . We also add the multiplicative 
shear corrections by replacing 7 with (1 -|- 771 ) 7 . We note 
that the observed ellipticity for a particular galaxy al¬ 
ready contains the leasing shear by large scale structure 
(LSS), but the random rotation makes this contribution 
at least second order in k by destroying the shape spa¬ 
tial correlations induced by leasing from LSS. Consistent 
with the weak leasing approximation, the leasing signal 
from the simulations is first order in k and hence the 
randomly rotated observed ellipticities can be safely con¬ 
sidered as intrinsic ellipticities. 

We analyze the simulations in the same way as we ana¬ 
lyzed the CFHTLenS data - constructing the simulated k 
maps as explained in llll Al These final simulation prod¬ 
ucts are then processed together with the k maps ob¬ 
tained from the data to compute confidence intervals on 
the parameter triplet (Hm, w,erg)- 


III. STATISTICAL METHODS 

The goal of this section is to describe the framework to 
combine the CFHT data and our simulations, and to de¬ 
rive the constraints on the cosmological parameter triplet 
{flm, w, CTs)- Briefly, we measure the same set of statisti¬ 
cal descriptors from the data and from the simulations; 
these are then compared in a Bayesian framework in or¬ 
der to compute parameter confidence intervals. 


A. Descriptors 

The statistical descriptors we consider in this work 
are the Minkowski Functionals (MFs) and the low-order 
moments (LMs) of the convergence field. The three 
MFs (Vq, Fi, V 2 ) are topological descriptors of the conver¬ 
gence field k{0), probing the area, perimeter and genus 
characteristic of the k excursion sets defined as 

Sko = {k > kq}- Following refs. [H, [ 2 ^ we use the 
following local estimators to measure the MFs from the 
K maps: 

K)(ko) = ^ y - Ko)M, 

Vi{ko) = + KldB, ( 8 ) 


^ 2 (^ 0 ) — 


27rA 


5d{k{0)-kq) 


2^Kxl^yl^Xy ^x^yy ^y^X 


Here A is the total area of the field of view and Kx,y de¬ 
notes gradients of the k field, which we evaluate using 
finite differences. In this notation 0(a;) is the Heaviside 
function and Sd{x) is the Dirac delta function. The first 
Minkowski functional, Vq, is equivalent to the cumulative 
one-point PDF of the k field, while Vi , V 2 are sensitive 
to the correlations between nearby pixels. The one-point 
PDF of the K field, i9Foi can be obtained by differentia¬ 
tion dVo{KQ) = dVo{tio)/dKo. 

In addition to these topological descriptors, we con¬ 
sider a set of low-order moments of the convergence field 
(two quadratic, three cubic and four quartic). We choose 
these moments to be the minimal set of LMs necessary to 
build a pert urbative expansion of the MFs up to O(o'g) 
(see [4l|,|^). We adopt the following definitions 

LM 2 :a2i = («2),(|V7c|2), 


LM3 : 5 o.1,2 = (k^), (k|Vk|2), (k2v2k). 


(9) 


LM4 : ifo.i.2.3 = (|Va«| 4 ). 

If the K field were Gaussian, one could express the MFs 
in terms of the LM 2 moments, which are the only inde¬ 
pendent moments for a Gaussian random field. In re¬ 
ality, weak lensing convergence fields are non-Gaussian 
and the MF and LM des crip tors are not guaranteed to 
be equivalent. Refs. studied a perturbative ex¬ 

pansion of the MFs in powers of the standard deviation 
(To of the K field. When truncated at order 0 ((Tq), this 
can be expressed completely in terms of the LMs up to 
quartic order. Such perturbative series, however, have 
been shown not to converge [ 2 ^ unless the weak lensing 
fields are smoothed with windows of size > 15'. Because 
of this, throughout this work, we treat MF and LM as 
separate statistical descriptors. 

We note that this choice is somewhat ad-hoc. In gen¬ 
eral, the LMs that contain gradients are sensitive to dif¬ 
ferent shapes of the k multispectra P"(li,..., 1 „) because 
a particular LM„ has the general form 

LM„ = J dh...dlnPill...n)P:{h...n) ( 10 ) 

where p is a polynomial of order n in the I’s. For exam¬ 
ple for K 2 we have p(li 234 ) = ^4 and this moment em¬ 
phasizes quadrilateral shapes for which one side is much 
larger than the others. On the other hand, for we 
have p(li 234 ) = (li • l 2 )(l 3 • I4) and this moment is most 
sensitive to trispectrum shapes that are close to rectan¬ 
gular. There are moments which include derivatives in 
addition to those included in Eq. O In the future, we 
will investigate whether there is additional constraining 
power in these additional quartic moments. 

In addition to the MFs and LMs, we consider the an¬ 
gular power spectrum Pi = of k, defined as 

(It(l)k(l')) = (2^)25^(l + l')Pi, (11) 


d0. 
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Descriptor 

Details 

Nb (linear spacing) 

Uo,Ui,U 2 (ME) 

Ko G [-0.04,0.12] 

50 

Power Spectrum (PS) 

1 G [300, 5000] 

50 

Moments (LM) 

- 

9 


TABLE II. Summary of the descriptors we used, together with 
the specifications and the number of bins Nb in each case. 


where k( 1) is the Fourier transform of the k field. Previ¬ 
ous works have studied cosmological constraints from the 
convergence power spectrum extensively. Here our pur¬ 
pose is to compare the constraints we obtain from the 
MFs and LMs to ones present in the literature, which 
are based on the use of quadratic statistics (see for ex¬ 
ample Q)- The statistical descriptors used in this work 
are summarized in Table im 

When measuring statistical descriptors on k maps, par¬ 
ticular attention must be paid to the effect of masked 
pixels. The MFs and LMs remain well-defined in the 
presence of masks, since the estimators in equations dH]) 
and dH]) are defined locally, and can be computed in the 
non-masked regions (with the exception of the few pix¬ 
els that are close to the mask boundaries). The situa¬ 
tion is more complicated for power spectrum measure¬ 
ments, which require the evaluations of Fourier trans¬ 
forms and hence rely on the value of every pixel in the 
map. Although sophisticated schemes to interpolate over 
the masked regions have been studied (see for example 
[ 4 ^ 1 . for the sake of simplicity, we here insert the value 
K = 0 in each masked pixel. Given the uniform spa¬ 
tial distribution of the masked regions in the data, we 
expect that masks have a little effect on the power spec¬ 
trum at the range of multipoles in Table m except for 
an overall normalization which will be the same both in 
the data and the simulations. Likewise, we believe that 
the way we deal with masked sky regions - essentially ig¬ 
noring them - is robust for the MFs and LMs. Since we 
apply the same masks to our simulations and the data, 
they are unlikely to introduce biases in the resulting con¬ 
straints. Masks, of course, can still affect the sensitivity 
and weaken constraints. The impact of the masks and 
their treatment has been evaluated for the MFs, obtained 
from CFHTLenS, by ref. [ 2 ^, in which the authors find 
that the masked regions are not a dominant source of 
systematic effects in the CFHTLenS data. 


B. Cosmological parameter inferences 

In this section, we briefly outline the statistical frame¬ 
work adopted for computing cosmological parameter con¬ 
fidence levels. We make use of the MFs and LMs, as well 
as the power spectrum, as discussed in the previous sec¬ 
tion. We refer to M[ (p) as the descriptor measured from 
a realization r of one of our simulations with a choice of 
cosmological parameters p (i.e. from one of the R = 1000 
map realizations in this cosmology), and to Di as the 


descriptor measured from the CFHTLenS data. In this 
notation, i is an index that refers to the particular bin on 
which the descriptor is evaluated (for example i can range 
from 0 to 9 for the LM statistic and from 0 to iVf, — 1 for 
a MF measured in Ni, different, linearly spaced k bins, 
as indicated in Table El). 

Once we make an assumption for the data likelihood 
J^d{Di\p) and for the parameter priors n(p), we can use 
Bayes’ theorem to compute the parameter likelihood Cp, 


_) ^ ^-(AMn(p) 


( 12 ) 


Here Nc is a p-independent constant that ensures the 
proper normalization for Cp. We make the usual assump¬ 
tion that the data likelihood Cd{Di\p) is Gaussian 


Cd{D,\p) = [(27r)^'>detC]-i/2e-^x"(o.|p)^ 

(13) 

X2(A|P) = [D - M(p)]"’c- 1[D - M(p)]. 

We assume, for simplicity, that the covariance matrix 
C in equation (1131) is p-independent and coincides with 
C(po). The simulated descriptors M(p) are measured 
from an average over the R = 1000 realizations in the 
CFHTemul ensemble 


M,(p) 


1 


R 

R^ * 

r-=l 


(14) 


The covariance matrix 
1 ^ 

C'b = Et^^Po) - M,(po)][MJ(po) - M,(po)] 

r—1 

(15) 

is measured from the R = 1000 realizations in the 
CFHTcov ensemble. While eq. dig gives an unbiased es¬ 
timator of the covariance matrix, its inverse is not an 
unbiased estimator of (e.g. ref. [l^). Given that in 
our case i? » A, we can safely neglect the correction 
factor needed to make the estimator for C~^ unbiased. 

When computing parameter constraints from the 
GFHTLenS weak leasing data alone, we make a flat prior 
assumption for n(p). We postpone using different priors, 
incorporating external data, for future work. Parameter 
inferences are made estimating the location of the max¬ 
imum of the parameter likelihood in eq. which we 
call PMbiDi), as well as its confidence contours. The 
A^cr-confidence contour of £p{p\Di) is defined to be the 
subset of points in parameter space on which the likeli¬ 
hood has a constant value cn and 

p 1 pN 

/ Cp{p\D,)dp = —j= / (16) 

JOcN V 27r J-N 

Using equation (13 below, and given the low dimen¬ 
sionality of the parameter space we consider [D = 3), 
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we are able to directly compute the parameter likeli¬ 
hood eq. m for 100^ different combinations of the 
cosmological parameters p, arranged in a finely spaced 
100 X 100 X 100 mesh within the prior window n(p). We 
directly compute the maximum likelihood puhiDi) and 
the contour levels cn without the need for more sophis¬ 
ticated MCMC methods. 

The data likelihood is directly available for parame¬ 
ter combinations on the simulated irregular grid Ps- We 
use a Radial Basis Function (RBF) scheme to interpolate 
M(p) to arbitrary intermediate points. We approximate 
the model descriptor as 

N 

^(p) = 

s=l 

where (j) has been chosen as a multi-quadric function 
(j){r) = y/l + (r/ro)^, with ro chosen as the mean Eu¬ 
clidean distance between the points in the simulated grid 
Ps. The constant coefficients As can be determined by 
imposing the N constraints M(p = Ps) = M(ps), which 
enforce exact results at the simulated points. The in¬ 
terpolation computations are conveniently performed us¬ 
ing the interpolate. Rbf routine contained in a Scipy 
library [45l| . 

We studied the accuracy of the emulator, built with 
the CFHTemul simulations, by interpolating the con¬ 
vergence descriptors to the fiducial parameter setting 
{nm,w,<Ts) = (0.26,—1.0,0.8) and comparing the result 
to the one expected from the CFHTcov simulations. Fig¬ 
ure [2] shows that our power spectrum emulator has a 
relative error smaller than 20% for the lower multipoles 
{I < 500), and comparable to 1% for the higher mul¬ 
tipoles. The MF emulator has a relative error <10% 
for the first 30 bins (which correspond to k values in 
[—0.04, 0.08]) and deteriorates due to numerical noise for 
the remaining 20 bins. We eliminate the residual im¬ 
pact of these inaccuracies using a dimensionality reduc¬ 
tion framework, which we explain in the next section. 
Nevertheless, we do not expect these inaccuracies to af¬ 
fect our conclusions. Figure [2] demonstrates that our em¬ 
ulator is able to distinguish a non-fiducial model from 
the fiducial one within numerical errors. We thus found 
no need to implement a more sophisticated interpolation 
scheme [i^ . 

C. Dimensionality reduction 

The main goal of this work is constraining the 
cosmological parameter triplet (flm, n;, erg) using the 
CFHTLenS data. Once the N<t contours have been ob¬ 
tained, using the procedure and equations dT^ - lfTHI) out¬ 
lined above, one may ask whether the choice of binning 
affects these contours. Indeed, in our previous work, we 
have found that the number of bins, W, can have a non- 
negligible effect on the contour sizes (see for an ex¬ 
ample with simulated datasets). 


In order to ensure that are results are robust with re¬ 
spect to binning choices, we have implemented a Prin¬ 
cipal Component Analysis (PCA) approach. Our physi¬ 
cal motivation for this approach is that, even though we 
need to specify Nt numbers in order to fully characterize 
a binned descriptor, we suspect that the majority of the 
constraining information (of a particular descriptor) is 
contained in a limited number of linear combinations of 
its binning. In the framework adopted by [3^ . for exam¬ 
ple, the authors find that the majority of the cosmological 
information in the matter power spectrum is contained 
in only 5 different linear combinations of the multipoles. 
Because of this, we believe that dimensionality reduc¬ 
tion techniques such as PCA can help deliver accurate 
cosmological constraints using only a limited number of 
descriptor degrees of freedom. 

In order to compute the principal components of our 
descriptor space, we use the CFHTemul simulations, which 
sample the cosmological parameters at the = 91 points 
listed in Table HI and allow us to compute the N x Nb 
model matrix Mpi = Mi{p). Note that this is a rectangu¬ 
lar (non-square) matrix. Following a standard procedure 
(see, e.g., ref. |^), we derive the whitened model matrix 
Mpi, defined by subtracting the mean (over the N = 91 
models) of each bin, and normalizing it by its variance 
(always over the iV = 91 models). Next we proceed with 
a singular value decomposition (SVD) of M, 

usv^ = , (18) 

where Sij = SiSij is a diagonal matrix and is the j-th 
coordinate (j = l...Nb) of the i-th principal component 
(i = I...min[Af,, N]) of M, with the index j ranging from 
I to Aft. By construction, V is p-independent. 

To rank the Principal Components in order of im¬ 
portance, we note that the diagonal matrix is simply 
the diagonalization of the model covariance (not to be 
confused with the descriptor covariance in eq. (1X5]) ]). 

M^M = VS^V^. (19) 

We follow the standard interpretation of PCA compo¬ 
nents, stating that the only meaningful components 
in the analysis (i.e. the ones that contain the relevant 
cosmological information) are those corresponding to the 
largest eigenvalues Sf, with the smallest eigenvalues cor¬ 
responding to noise in the model, due to numerical inac¬ 
curacies in the simulation pipeline. We expect our con¬ 
straints to be stable with respect to the number of com¬ 
ponents, once a sufficient number of components have 
been included. Using the fact that different principal 
components are orthogonal, we perform a PCA projec¬ 
tion on our descriptor space by whitening the descriptors 
and computing the dot product with the principal com¬ 
ponents, keeping only the first n components 

M(n)F = U^(n)yMF ; D{n), = {n),,b,. (20) 
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FIG. 2. Accuracy of the emulator based on the CFHTemul simulations. The figure shows the absolute difference between the 
descriptor interpolated at the fiducial parameter setting, and the descriptor expected from the CFHTcov simulations (these are 
the absoulte values of differences which oscillate around zero). The descriptors are shown in units of the standard deviation in 
each bin i (determined from the diagonal elements of the CFHTcov covariance matrix). We show the accuracy results for the 
power spectrum (red) and the three Minkowski functionals Vo (green), Vi (blue) and V 2 (black). For reference, we also show, 
using dashed lines, the difference between the expected CFHTcov descriptors and the interpolated descriptor at the non-fiducial 
point p = (0.8, —1.0, 0.5). This non-fiducial point lies beyond the LM Icr contour from the simulations shown in Figure |3] right 
panel, and corresponds to the target accuracy we wish to achieve 


Here we indicate with V^(n) the truncation of V'^ to 
the first n rows (i.e. i can now range from 1 to n). As 
described above, the expectation is that most of the cos¬ 
mological information is contained in a small number of 
components n < Nb- We will describe in detail below the 
choice we make for n, together with the sensitivity of our 
results to this choice. 

Looking at PCA from a geometrical perspective, the 
dimensionality reduction problem is equivalent to the ac¬ 
curate reconstruction of the coordinate chart of the de¬ 
scriptor manifold. As outlined in ref. [d^, the coordinate 
chart constructed with the PCA projection in eq. (1201) is 
accurate for reasonably flat descriptor manifolds. When 
curvature becomes important, more advanced projection 
techniques (such as Locally Linear Embedding) have to 
be employed. We postpone an investigation of such im¬ 
provements to future work. 

IV. RESULTS 

This section describes our main results, and is orga¬ 
nized as follows. We begin by showing the cosmologi¬ 
cal constraints from the CFHTLenS data for the triplet 
(Hm,<78, w), as well as for an alternative parameteriza¬ 
tion, (Hm,E8,n'), with ^ 8 ( 0 :) = (78(^^m/0.27)“. In the 
next section we give a justification on why we fix a to a 
value of 0.55. We then use our simulations to perform a 


robustness analysis of the parameter confidence intervals 
with respect to the number of PCA components used in 
the projection. We finally study whether the constraints 
can be tightened by combining different descriptors. A 
summary with the complete set of results, along with the 
relevant Figures, is shown in Table Hill 

A. Cosmological constraints 

We first make use of equations (fT2l) - (|T6)) to compute 
the la constraints on cosmological parameters, using 
the triplet {Qm,as,'w). Figures [3] shows the constraints 
in the (flrmcra) plane, marginalized over w, for both 
the CFHTLenS data, as well as from the mock data in 
our simulations. In Figure HI we examine constraints 
from different sets of moments, as well as using different 
smoothing scales. Figure [S] shows the confidence contours 
in the (w, E8) plane, marginalized over 17^. As this fig¬ 
ure shows, and as discussed further below, no meaningful 
constraints were found on w from CFHTLenS alone. 

Because of the relatively small size of this survey, de¬ 
generacies among the parameters can have undesirable 
effects on the constraints. The well-known strong de¬ 
generacy between Qm and as is evident in the long “ba¬ 
nana” shaped contours in Figures [3] and HI To mitigate 
the effect of this degeneracy, in addition to the usual 
triplet (Hm, 178, re), we consider an alternative parame- 
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terization, built with the triplet (a, Sgji/;) where a is a 
constant, and ^s{c() = o'8(nm/0-27)“. While and 
tTs are poorly constrained due to degeneracies, the S8(a) 
combination lies in the direction perpendicular to the er¬ 
ror “banana” at the pivot point ilm = 0.27. This is the 
direction of the lowest variance for a suitable 

choice of a, and hence has a much smaller relative un¬ 
certainty. We can derive the optimal value of a from 
the full three dimensional likelihood C{Qrn,WT crs), from 
which we can compute the expectation values 

E(a) = (Esia)) ; V(a) = ((EgCa) - E(a))2) (21) 

and minimizing the ratio /E with respect to a. The 
expectation values are taken over the entire parameter 
box. This procedure yields a value a ~ 0.55 for the sta¬ 
tistical descriptors that we consider, consistent with what 
is found in the literature (see [ij for example). Although 
a can mildly depend on the type of descriptor consid¬ 
ered, we choose to keep it fixed, knowing that the width 
of the S8 likelihood cannot vary significantly with differ¬ 
ent choices of a. We show the probability distribution 
of the best-constrained parameter Es (marginalized over 
flm and w) in Figure [H 

We discuss the results of this section in § |V] below. 


B. Robustness 

The cosmological constraints should in principle be in¬ 
sensitive to A^b, once a sufficient number of bins are used, 
but inaccuracies in the covariance (due to a limited num¬ 
ber of realizations) can introduce an N\y dependence. Our 
binning choices are summarized in Tableln] Here we show 
that the cosmological constraints derived in this paper 
are numerically robust, i.e. they are reasonably stable, 
once we consider a large enough number n of Principal 
Components. 

Figure [7] shows the PC A eigenvalues from the SVD 
decomposition of our binned descriptor spaces (follow¬ 
ing the discussion in § IIII Cp , as well as the cumulative 
sum of these eigenvalues, normalized to unity. Figure [5] 
shows the dependence of the (Hm, crs) constraints on the 
number of principal components n. 

These figures clearly indicate that we only need a lim¬ 
ited number of components in order to capture the cos¬ 
mological information contained in our descriptors. The 
eigenvalues diminish rapidly with n, and, in particular, 
the confidence contours converge to good (< 10%) accu¬ 
racy typically for n = 5 — 10 (depending on the descrip¬ 
tor). This finding also addresses the inaccuracy of the 
MF emulator at high thresholds, pointed out in Figured 
By keeping a limited number of principal components, we 
are able to prevent the inaccurate high-threshold bins, 
which have a low constraining power, from contributing 
to the parameter confidence levels. 


C. Combining statistics 

Different statistics can include complementary cos¬ 
mological information, allowing their combinations to 
tighten the constraints. Previous work using multiple 
lensing descriptors in CFHTLenS alone included combin¬ 
ing the power spectrum and peak counts [2^ , combining 
the power spectrum and Minkowski functionals ,and 
combining quadratic (2PCF) statistics with cubic statis¬ 
tics derived from the 3PCF of the CFHTLenS k field [^ . 

The procedure we adopt here is as follows. Con¬ 
sider two binned descriptors, di,i,d 2 j where the indices 
correspond to bin numbers. We first compute each 
single-descriptor constraint as a function of the number 
of PC A components, as in Figure [5] We then determine 
the minimum number of PCA components ni _2 needed 
for the constraints to be stable. We next construct the 
vector dix 2 = {di(rii),d 2 {n 2 )} and consider this as the 
combined (ni -|- n 2 )“dimensional descriptor vector. This 
procedure naturally allows us to account for the cross¬ 
covariance between different binned descriptors. An anal¬ 
ogous procedure can be used to combine multiple (three 
or more) descriptors. 

We show constraints from different descriptor combi¬ 
nations in the (Dm,o'8) and (re, Es) planes in Figure IHl 
and on the best-constrained parameter Es in Figure [TUI 
We also provide a tabulated version of the Es constraints 
(Icr) in Table llVl We discuss these findings in the next 
section. 


V. DISCUSSION 

In this section we discuss the results shown in §IIY] 
above, with particular focus on the constraints on cos¬ 
mology. 

As pointed out in § IIII Cl the choice of the number 
of bins, Nfn is an important issue. In order to ensure 
that our results are insensitive to Nt, we adopted a PCA 
projection technique to reduce the dimensionality of our 
descriptor spaces. The left panel of Figure |7| shows that 
the PCA eigenvalues for all of our descriptors decrease by 
about 4 orders of magnitude from n = 1 to n = 3. The 
right panel of this figure shows that more than 99% of 
the descriptor variances are captured by including only 
the first n = 3 components. 

This does not necessarily mean, however, that the cos¬ 
mological information is captured by the first 3 PCA 
components: in principle, one of the higher-n PCA 
components could have an unusually strong cosmology- 
dependence, and could impact the confidence levels. To 
address this possibility, we determined the Icr contour 
sizes as a function of n in Figure This figure shows 
that the first 3 components indeed capture essentially all 
the information contained in the power spectrum. How¬ 
ever, this is not true for the other descriptors. In partic¬ 
ular, we find that n > 5 components are necessary for Vq, 
and n > 20 components for Vi and V 2 , in order for the 
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FIG. 3. la (68% CL) constraints on the (fImjCrg) parameter doublet using the power spectrum (red), the three Minkowski 
functionals (Vo: green, Vi: blue, V 2 : black) and the moments (orange). We show the constraints from the data (left panel) and 
from a mock observation constructed using the mean of 1000 realizations in the CFHTcov simulation suite (right panel). The 
contours are calculated from the parameter likelihood function £ marginalized over w. The parentheses near the descriptor 
label refer to the number of principal components included. 




FIG. 4. Icr (68% CL) constraints on the {flm,as) parameter doublet using moments, with different colors corresponding to 
different moment combinations (see eq. [9]for their definitions). We show the results from the one-point moments aQ,So,Ko 
(black curves; both left and right panels). In the left panel, we also show constraints obtained adding moments of gradients to 
the one-point moments. In the right panel, we combine one-point moments measured at different smoothing scales. 
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Parameters 

Descriptors 

Short description 

Relevant Figures 

(r^TTi, (Tg) 

PS(3),yo(5),Fi(20),F2(20),LM(9) 

Icr constraints from GFHl'LenS 

and mock observations 

[3l3b 

(r^TTi, fjg) 


Icr constraints from Cf HTLenS 

using K moments 
combined at different ffn 

I4l4b 

(ly. Eg) 

PS(3),V'o(10), Vi(lO), V2(10),LM(9) 

contours from CFHTLenS 

0 

Es 

PS(3),Vo(10), Vi(lO), V2(10),LM(9) 

£(Eg) from CFHTLenS 

0 

- 

PS,Vo,Fi,F2,LM 

PCA eigenvalues 

0 

(r^TTi, fjg) 

PS,Vo,Vi,V2,LM 

Stability of contours 

0 

(r^TTi, fTg) 

PS(3)xVo(5) X Vi(20) X V2(20)xLM(9) 

constraints from CFH'i'LenS 

combining statistics 

0 

(ry.Eg) 

PS(3)xVo(10) X Vi(lO) X F2(10)xLM(9) 

constraints from CFH'i'LenS 

combining statistics 

0D 

Eg 

PS(3)xVo(10) X Vi(lO) X V2(10)xLM(9) 

£(Eg) from CFH'i’LenS 
combining statistics 

m 


TABLE III. Summary of our results and related figures. 


Descriptor 

Eg = agH^®® 

PS(3) 

PS(3) X Moments(9) 

Vo(lO) X Vi(lO) X V2(10) 

PS(3) xVo(lO) X Vi(lO) X V2(10) 

PS (3) xVo(lO) X Fi(lO) X V2(10)x Moments(9) 

0.84t°;°® 

0.861^9 

0.75l°;°I 

n 7^+0 04 

U. ^O_o.o5 

0 76+0 06 


TABLE IV. Tabulated values of la constraints on Es corresponding to Figure [TO] 



FIG. 5. la (68% C.L.) constraints on the (WjEs) parameter 
doublet from the CFHTLenS data, obtained with the power 
spectrum (red), the three Minkowski functionals (Vq: green, 
Vi: blue, Vi: black) and the moments (orange). The con¬ 
tours are calculated from the parameter likelihood function £ 
marginalized over Q^n, and the parentheses near the descrip¬ 
tor label refer to the number of principal components. 


{flm, o's) contours to be stable at the ~5% level. All nine 
moments need to be included for the moments contours 
to be stable to this accuracy. 

These results are slightly different when we study the 
(w. Eg) constraints (with a fixed at a = 0.55 as discussed 
above). In this case we find that the optimal choice for 
all three MFs is n = 10, while the number of components 
required for the PS and Moments remain at n = 3 and 
n = 9, respectively. (These results are not shown, but 
obtained analogously to the Figures above.) 

We now discuss the main scientific findings of this 
work. In Figure [Sj we show the Icr constraints on 
the (firnjO's) doublet from the CFHTLenS data. The 
MF constraints appear to be biased towards the low- 
CTg, high-Om region. Here and throughout the remain¬ 
der of this paper, by ’’biased” (or ’’unbiased”) we re¬ 
fer to being incompatible (or compatible) with the con¬ 
cordance fiducial values at Icr obtained in other ex¬ 
periments. For example, the current best-fit values 
of (rimjO-g) = (0.32, 0.83), (0.28,0.82) from cosmic mi¬ 
crowave background anisotropies measured respectively 
by the Planck and WMAP [s^ satellites lie beyond 
the 99% likelihood contours obtained from the three MFs 
(not shown). 

This discrepancy may be due to uncorrected system- 
atics in the CFHTLenS data, amplified by the (HmjCrg) 
degeneracy. As a test of our analysis pipeline, when we 
try to constrain mock observations based on simulations 
(shown in the right panel of Figure[3|), we recover the cor¬ 
rect input position of the Icr contours. It is important to 
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FIG. 6. The likelihood of the best-constrained parameter 
combination E8 (q;) = (T8(flm/0.27)“ from the CFHTLenS 
data, obtained with the power spectrum (red), the three 
Minkowski functionals (Vq: green, Vi: blue, V 2 : black) and 
the moments (orange). The likelihood was computed with a 
constant optimized a = 0.55, but marginalized over both flm 
and w. The parentheses near the descriptor label refer to the 
number of principal components. 


note, however, that the mock observations to which the 
right panel of Figure [3] refers, were built with the mean 
of i? = 1000 realizations of the CFHTcov simulations. We 
found that it is possible to find some rare (10 out of the 
1 ,000) realizations for which the best fit for (flmjfs) lies 
in the lower right corner, near the location of the best-fit 
from the data. While this could provide an alternative 
explanation of the bias from the MFs, the likelihood of 
this happening is very small (< 1%). 

We observe that the moments give the tightest con¬ 
straint on (flmjCrg). Furthermore, this constraint is un¬ 
biased, in the sense defined above: it includes the current 
concordance values for these parameters within Icr. This 
leads us to conclude that the bias in the constraints from 
the MFs is due to systematic errors, rather than the rare 
statistical fluctuations found above. The fact that the 
moments are useful for deriving unbiased cosmological 
constraints has been noted in previous work, which ex¬ 
amined the biases caused by spurious shear errors [i^ . 

In order to determine the origin of the tight bounds de¬ 
rived from moments, we studied the contribution of each 
individual moment to the constraints. Figure S] shows the 
evolution of the (flmiCrs) constraints as we add increas¬ 
ingly higher-order moments to the descriptor set. Since 
we are constraining 3 cosmological parameters, we start 
by considering the set of the three traditional one-point 


moments which do not involve gradients, i.e. the vari¬ 
ance, skewness, and kurtosis (aQ, So, Kg). We then add 
the remaining six moments of derivatives one by one, 
starting from the quadratic moments. 

Figure 0] shows that the biggest improvement on the 
parameter bounds comes from including quartic moments 
of derivatives (i.e. Ki with i > 1) in the descriptor set. 
This might explain why find only relatively weak 
contour tightening (~ 10%) when adding three-point 
correlations to quadratic statistics, since the main im¬ 
provement comes from higher moments of k derivatives. 
Ref. consider one-point, third-order moments, com¬ 
bined for multiple smoothing scales. Figure 0] explicitly 
shows, however, that smoothing scale combinations are 
not as effective as moments of derivatives in constraining 
the doublet. Our results agree with an early 

prediction Q that the kurtosis of the shear field can help 
in breaking degeneracies between flm and erg. Here we 
found that considering quartic moments of gradients fur¬ 
ther helps in breaking this degeneracy. 

As noted above, the bias in the (0^, erg) constraints is 
amplified by the cosmological degeneracy of these param¬ 
eters. To mitigate this effect, we consider the combina¬ 
tion of fim and tTg that lies orthogonal to the most degen¬ 
erate direction, namely Sg = crg(Om/0.27)°'®^. Figure^ 
shows the Icr constraints for the (w, Eg) doublet, while 
Figure [6] shows the marginalized Eg likelihood from the 
CFHTLenS data. The CFHTLenS survey constrains the 
Eg combination to a value of Eg = 0.75 ± 0.04(lcr) using 
the full descriptor set, in agreement with values previ¬ 
ously published by the CFHTLenS collaboration [l|. 

These figures also show that the current dataset is in¬ 
sufficient to constrain w to a reasonable precision. This 
is consistent with the previous analyses of CFHTLenS [l|, 
m, HE HI] • We also note that ref. obtained the best- 
fit value of w ~ —2 (but with large errors that include 
w = —1 at Icr). We found a similar result when using 
a Fisher matrix to compute confidence levels. Since the 
Fisher matrix formalism is equivalent to a linear approx¬ 
imation of our emulator (in which all cosmological pa¬ 
rameter dependencies are assumed to be linear), we thus 
attribute this bias to the oversimplifying assumption of 
linear cosmology-parameter dependence of the descrip¬ 
tors. Although the right panel in Figure [9] shows that 
the moments confine w to isolated regions in parameter 
space, we note that w = —1, the value favored by other 
existing experiments, is excluded at the Icr level. The 
2 cr contours (not shown in the figure) join, and include 
w = —1. 

Regarding the parameter biases, our results overall 
are in accordance with (^ . namely, that unaccounted 
systematics result in larger parameter biases when the 
constraints are derived from the MFs, and that the LM 
statistic is less biased. However, for the CFHTLenS data 
the MFs can still effectively constrain the non-degenerate 
direction in parameter space, Eg (Figure 0]). 

Finally, we studied whether the combination of dif¬ 
ferent statistical descriptors can help in tightening the 
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FIG. 7. Results from the Principal Components Analysis (PCA) of the binned power spectrum (red), the three Minkowski 
functionals {Vq: green, Vi: blue, Vi: black) and the moments (orange). The left panel shows the magnitudes of the PCA 
eigenvalues Sf and the right panel shows their cumulative sum, normalized to unity. A dashed vertical (black) line has marks 
n = 3 for reference. 


cosmological constraints. We show the effects of some of 
these combinations in Figures l^andfTUl The left panel of 
Figure |9] shows that, although combining the power spec¬ 
trum and the moments with the Minkowski functionals 
helps tighten the (flmjCrg) constraints, it does not help 
in reducing the inherent parameter bias of the MFs. The 
right panel of Figured shows that even with these statis¬ 
tics combined, w remains essentially unconstrained, 
found that even weak leasing tomography alone is un¬ 
able to constrain w sensibly. Figure [in] shows that the 
Sg combination is already well constrained by any of the 
descriptors alone, without the need of combining different 
descriptors. This further clarifies that the non-quadratic 
descriptors mainly help to break degeneracies, tightening 
contours along the degenerate direction. 


VI. CONCLUSIONS 

In this final section we summarize the main conclusions 
of this work: 

• We find that the power spectrum, combined with 
the moments of the n held provides the tight¬ 
est constraint on the doublet from the 

CFHTLenS survey data. The tightness of these 
constraints comes mainly from the moments. Evi¬ 
dence of the unbiased nature of constraints from the 
moments has been found in [d^. We further hnd 
that the largest improvement on parameter bounds 
is achieved when we include the quartic moments of 
derivatives in the descriptor set. This level of im¬ 


provement cannot be achieved by combining one- 
point moments at different smoothing scales. 

• Although weak leasing surveys are a promising 
technique to constrain the DE equation of state 
parameter w, reasonable constraints cannot be ob¬ 
tained with the CFHTLenS survey alone, even 
when using additional sets of descriptors that go 
beyond the standard quadratic statistics. 

• When studying the cosmological information con¬ 
tained in the CFHTLenS data, special attention 
must be paid to the effect of residual systematic 
biases. While these residual systematics are found 
to be unimportant when constraining cosmology 
with the power spectrum alone, we find that these 
systematics need to be corrected to obtain unbi¬ 
ased constraints on the (Hm,CTg) doublet using the 
Minkowski functionals. We are aware that, when 
trying to explain the discrepancy between weak 
leasing and CMB constraints using the Minkowski 
functionals, there might be other effects to be con¬ 
sidered, namely non-Gaussian error correlations in 
the descriptors and inaccuracies of the simulations 
on small scales. These inaccuracies could in princi¬ 
ple affect the excursion set reconstruction at high 
kq thresholds. We will investigate these additional 
sources of error in future work. 

• For the CFHTLenS data set, Minkowski functionals 
can effectively constrain the non-degenerate direc¬ 
tion in parameter space. Eg, where the amplifying 
effects of degeneracy are mitigated.The Minkowski 
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FIG. 8. The dependence of the Icr contours in the (flmjCrs) plane on the number of PCA components, obtained from a mock 
observation constructed with the CFHTcov simulations. The different panels refer to the different descriptors (from left to right, 
top to bottom) Vo, 9Vb(PDF), Vi, Vo, power spectrum and moments. The labels in each panel show the number of PCA 
components included to obtain contour with different colors. 
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FIG. 9. 1(7 constraints on the (fl^jCrg) (left panel) and (w, Eg) (right panel) doublets, using the power spectrum (PS) alone 
(red), the MFs alone (blue), as well as using different combinations of descriptors: PSxMoments (green), PSxMFs (black) and 
PSxMFsxMoments (orange). The likelihood function has been marginalized over w (left panel) and Qm (right panel). The 
parentheses next to each descriptor label refers to the number of PCA components included. 



FIG. 10. The probability distribution of the best-constrained 
parameter Eg from the GFHTLenS data, using the power 
spectrum (PS) alone (red), the MFs alone (blue), as well 
as using different combinations of descriptors: PS x Moments 
(green), PSxMFs (black) and PSxMFsxMoments (orange). 
The likelihood function has been marginalized over Qrn and 
w. The parentheses next to each descriptor label refers to the 
number of PCA components. 


functionals alone are sufficient to constrain the Eg 
combination to a value of Eg = 0.75 ± 0.04 at Icr 
significance level. This agrees with the value pre¬ 
viously published by the GFHTLenS collaboration 
within IfT. Some tensions with Planck still re¬ 
main. 

Possible future extensions of this work include simu¬ 
lating higher-dimensional parameter spaces (including, 
for example, the Hubble constant Hq, and allowing a 
tilt in the power spectrum or a time-dependence of the 
DE equation of state w), and combining the GFHTLenS 
constraints with different cosmological probes from large- 
scales structures and the GMB. The latter can help in 
breaking the (7g degeneracy, and allow improvements 
in the the constraints on w. The techniques developed 
here can be applied to larger, soon-forthcoming su rvey 
data sets, such as the Dark Energy Survey (DES) [5ll |. 
Subaru [Hi, WFIRST [H and LSST [5l|. 
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